Relaxation of hole spins in quantum dots via two-phonon processes 
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We investigate theoretically spin relaxation in heavy hole quantum dots in low external magnetic 
fields. We demonstrate that two-phonon processes and spin-orbit interaction are experimentally 
relevant and provide an explanation for the recently observed saturation of the spin relaxation rate 
in heavy hole quantum dots with vanishing magnetic fields. We propose further experiments to 
identify the relevant spin relaxation mechanisms in low magnetic fields. 
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In the last decade, remarkable progress has been made 
in the manipulation and control of the spin of electrons 
confined in semiconducting nanostructures such as quan- 
tum dots (QD) Q. These achievements pave the way 
toward quantum spin electronics and may lead to spin- 
based quantum computing j^,]. In the past years, a new 
candidate for a qubit state has been attracting growing 
interest: the spin of a heavy hole (HH) confined in a flat 
QD. In a bulk semiconductor the HH {J^ = ±3/2) and 
light hole (LH) (J^ = ±1/2) bands are degenerate giving 
rise to strong mixing and thus to strong HH-spin relax- 
ation. However, in a 2D system the HH and LH bands 
are split due to the strong confinement along the growth 
direction [3] implying a significant reduction of the HH 
spin relaxation via HH-LH mixing. 

The relaxation (Ti) and decoherence (T2) times of a 
HH-spin localized in a QD are, like for electrons, deter- 
mined by the environment the hole spin interacts with: 
the nuclear spin bath in the QD and the lattice vibra- 
tions (phonons). The former interaction is weaker for 
HHs than for electrons (due to the p-symmetry of the 
hole) ill. More importantly, it is of Ising type, mak- 
ing it ineffective for HH-spins initialized along the growth 
direction as typically done in experiments [61]. This 
is in contrast to electrons, where the hyperfine interac- 
tion is isotropic and dominates the spin dynamics at low 
B-fields 0,[ii,[ll. 

The other relevant source of relaxation are phonons 
which couple to the hole spin through the spin-orbit in- 
teraction (SOI) The predicted values pT| for the 
one-phonon induced relaxation time Ti agree quite well 
with data obtained in high B-fields . However, for low 
B-fields {B ~ 1.5-3T) and high temperatures (T > 2K), 
a clear deviation from the one-phonon theory has been 
observed Furthermore, recent experiments on opti- 
cal pumping of HH-spins in QDs showed saturation of Ti 
for very low or even vanishing B-field [6] . The relaxation 
time was found to be unusually long, Ti w 0.1 — 1ms, 
like previously observed in high B-fields [T^]. Both ob- 
servations suggest other sources of relaxation, and the 
question arises what are they and what are their observ- 
able consequences? The answer to this question is not 



only interesting by itself but also relevant for using HHs 
as qubits. In the following, we show that two-phonon 
processes are good candidates and even provide a quan- 
titative explanation of the mentioned measurements at 
low B-fields [T^]. Indeed, as we will see, these pro- 
cesses show weak or no B-field dependence, whereas one- 
phonon relaxation rates vanish quickly with B 0. The 
importance of such two-phonon processes was noticed a 
long time ago for electron spins in silicon-donors [l3| and 
rare-earth ions 14], while for electrons in QDs, they have 
been recently analyzed and shown to be negligible com- 
pared to nuclear spin effects J^]. 

To describe a HH confined to a QD and interacting 
with the surrounding phonon bath, we start with the 
following Hamiltonian 

Hh = Hq + Hz + Hso + Hh-ph + Hphi (1) 

where i/o = /Im* + V{r), is the dot Hamiltonian, 
V{r) = m*uo'^r'^ /2 is the confinement potential which 
is assumed to be harmonic, with m* being the HH mass. 
The second term in Eq. ([1]) is the Zeeman energy of 
the HH (pseudo-) spin Hz = gf^sB ■ cr/2, with B be- 
ing the magnetic field and cr the Pauli matrices for the 
HH spin defined in the Jz — ±3/2 subspace. The third 
term represents the spin-orbit Hamiltonian, which, for 
well separated HH-LH bands (flat dots), reads [ll[ 



Hso = PP-P+P-cr^ 



(2) 



This Hamiltonian represents the effective Dresselhaus 
SOI (restricted to the HH subspace) due to bulk inver- 
sion asymmetry of the crystal [ll| , where p± = Px± ipy , 
p = -ihV - eA{r), A{r) = {-y,x,0)B/2, and a± = 
ax ± ia-y We note that in Eq. ([2]) we have neglected the 
Rashba SOI and other possibly linear-in-k but small SOI 
terms fv\\. The fourth term in Eq. ^ represents the 
interaction of the HH charge with the phonon field, i.e. 



-ph = T,ni MqjXqj with 



F{qz 



i(^o Q ■ dqj - ^z Qzdq 



(3) 



and Xqj — yjh/ujqj{a}_q^ + flqj), whcrc q is the phonon 
wave-vector, with j denoting the acoustic branch, LOqj = 



2 



Cjq the phonon energy, with Cj the speed of sound in the 
j-th branch, dqj the polarization unit vector, pc the sam- 
ple density (per unit volume), and ePqj the piezoelectric 
electron-phonon coup ling and Sq ^ the deformation po- 
tential constants [ll|. The form factor F{qz) in Eq. ^ 
equals unity for <C and zero for \qz \ , with 

d being the dot size in the (transverse) z-direction. The 
last term in Eq. ([T]) describes the free phonon bath. 

In the following, we analyze the effect of the phonons 
on the HH spin, both in the low and high temperature 
regimes. The phonons do not couple directly to the spin, 
but the SOI plays the role of the mediator of an effective 
spin-phonon interaction. Under the realistic assumption 
that the level splitting in the dot is much larger than the 
HH-phonon interaction, we can treat Hh-ph in perturba- 
tion theory. 

Let us define the dot Hamiltonian Hd = Hq+Hz+Hso 
and denote by I^ct), the eigenstates of Hd (where n labels 
the orbital and a the spin states). The product states 
|n) |ct) are then the eigenstates of Hd in the absence of SOI 
{i.e. for Hso = 0). These states are formally connected 
to Incr) by an exact Schrieffer- Wolff (SW) transformation 
MEa. i-e-, \n<7) = e^\n)\a), where S = -S't is the SW 
generator and can be found in perturbation theory in 
SOI. After this transformation, any operator A in the 
old basis transforms as A ^ A = Ae~^ in the new 
basis (e.g., Hd Hd, Hh^ph Hh-ph, etc.). 

Let us now derive the effective spin-phonon interaction 
under the above assumptions. To do so, we perform an- 
other SW transformation of the total HH Hamiltonian 
Hfi to obtain an effective Hamiltonian H^jj — Hhe~^ , 
where T = — is chosen such that it diagonalizes the 
HH-phonon Hamiltonian Hh-ph in the eigenbasis of Hd- 
In lowest order in Hh-ph, we obtain T « L^^ Hh-ph, 
where the Liouvillean is defined as LdA = [Hd,A\, VA, 
and diagonal terms of Hh-ph are to be excluded. In 2nd 
order in Hh-ph, we obtain then the effective spin-phonon 
Hamiltonian 



^qlq'j' (^"Jj ^q'3' Pq'j' ^qj)] ' 



(4) 



with ^.C^yHO|M,,|0)^,T.C^^.;^,^HO|[i,-'M,,, M,,,.]|0), 

■ ^qlq'y = mLd'Mq„L^'Mq,,,]\0), Pq, = 

iy/^u^Jqjio^-qj — CLqj) IS the phouou field momentum 
operator, and |0) is the orbital ground state. In Eq. @ 
we have neglected spin-orbit corrections to the energy 
levels, being 2nd order in SOI. Note that for vanishing 
magnetic field B ^ the only term which survives 
in Hs-ph is the last one since this is the only one 
which preserves time-reversal invariance and thus gives 
rise to zero field relaxation (ZFR) [13, El, IH]. Quite 
remarkably, this term is 1st order in SOI, whereas for 
electrons it is only 2nd order [l5|. This is one of the 



main reasons why eventually two-phonon processes are 
much more effective for HHs than for electrons. 

We now assume the orbital confinement energy Uujq 
much larger than the SOI, i.e. \\Hq\\ ^ \\Hso\\, so that 
we can treat the SOI in perturbation theory. We consider 
also the B-field to be app lied perpendicularly to the dot 
plane (as in Refs. 0, [3). We limit our description to 
first order effects in SOI. The SW-generator S can be 
written as 5* = S+u- — h.c, and we then find 

5+= Aip+p-p+ + A2\p+p-P+ - [p+P- - P+p-)p+\ 
+ AiP+P-P+ + A^[{p+P- - P+p-)P+ + P+P-p+]{b) 

Here, Ai = Ai{ujz,^c) with loz ~ gfJ-BB/h and 
LUc ~ eB/2c. For luz,'-^c <C ujq, we obtain Ai « 
-i7l3/9n){{LOz A2 « ~il3/3h){uj,/uj^), A3 « 

^i2p/m){iu^^iu;, + iUz)/L0^), and A4 « i2f3/m)iu;^Ju;l), 
while P± ~ Px ± iPy with Px{y) — x(y) i 

{m*Lj'^/ijJc)y{x). After somewhat tedious calculations, we 
obtain analytic expressions for C^'^ = (C''^'^^ C^*'^', 0) 
occurring in Eq. (j4|). We give below only the exact 
expression for z = 3, the rest being too lengthy to be 
displayed here: 



C. 



qj,q ] 



- h^bx ± [hx - b'x){2hy% + ibX)) ,(6) 



where 



X (7 + log(6.672) + r(0,6-672)). (7) 



Here, M^/ = 



{F{q^)F{g',)h/2pcy/ujqjUjq,j>){Ea q ■ dq^j - 
qzdqjYi^o q' ■ dq'f - q'^d^q,j,) and b = qXd, where 
is the dot-diameter. We have also introduced 7 w 2.17 
the Euler constant and T{s,x) the incomplete gamma 
function. We note that C*^^'^-* cx B, so that these two 
terms vanish with vanishing B-field. 

Let us now analyze the relaxation of the spin induced 
by all the phonon processes in the spin-phonon Hamilto- 
nian in Eq. ([4]). We first mention that all terms in Eq. 
(HI can be cast in a general spin-boson type of Hamil- 
tonian = {l/2)gfiBSBP{t) ■ cT,p^ 1,2,3, with the 
corresponding identification of the fluctuating magnetic 
field terms 5Bj{t) from Eq. @ (e.g. 5B^{t) ~ &^^Xqj). 

Within the Bloch-Redfield approach, the relax- 
ation rate F = 1/Ti can be expressed as F = 
Y.i=x,y [J'liiEz/fi) + Jii{-Ez/h)]. The correlation func- 
tions Jij are defined by Jij{uj) — {g^s/'^K)^ /o°°^^s"^*< 
SBi{0)6Bj{t) >, where < • • • > denotes the average over 
the phonon bath, assumed to be in thermal equilibrium 
at temperature T. The relaxation time associated with 
the three types of spin-phonon processes in Eq. ^ is 



r(i) 

p(m) 



i=l,2,3 



r(') with 



An ■ 



(8) 



where ?T.(w) = l/(exp (w/fcsT) — 1) is the Bose factor and 
m — 2,3 correspond to B-dependent and B-independent 
two-phonon rates, resp. We remark that in Eq. ([8]) we 
have neglected some irrelevant processes in the limit of 
low-B field [l^]- Also, for B- fields perpendicular to the 
dot plane the decoherence time satisfies T2 = 2Ti for one- 
and two-phonon processes since the spin-phonon fluctu- 
ations (5Bj ± B 0,[T3. 

Note that for two-phonon processes the single phonon- 
energies do not need to match the Zeeman energy sep- 
arately (as opposed to one-phonon processes), so that 
there is only a weak dependence on the B-field left which 
comes from the effective spin-phonon coupling itself. 

In Figs. [1] and [21 we plot the phonon spin-relaxation 
rate F as a function of the B-field and of temperature, 
resp., for InAs and GaAs quantum dots. Fig. [1] shows 
a clear saturation of F at low magnetic fields which is 
due to two-phonon processes, while Fig. [2] shows the 
known saturation at low temperatures due to one-phonon 
processes 



m. 



For these plots, we used the following HH InAs QDs 
(labeled by A) [11, [2^ and GaAs QDs (labeled by B) 



parameters [11 
2.64 • 10^ m/s (c, 
(cf = 4.73 • 10 
5.3- lO^kg/m^), m 



1.9 eV, = 



2.7eV, cf = 



"0 

= 3.35 • 10^ m/s), = 3.83 • 10^ m/s 
m/s), = 5.68 • lO^kg/m^ (pf = 
= 0.25TOe {m*g ^ 0.14me), = 1.4 



{g^ = 2.5), and we assume — 3nm {Thuq — 35meV, 
Tiloq — 60meV) and d = 3nm (dot height). Also, 
Pa ~ 2.1 • 10^ m/s and /3b w 4.6 • 10"* m/s. From Fig. 
[1] we can infer that the two-phonon processes become 
dominant for magnetic fields B < 2T {B < 0.5T) and 
for temperatures T > 2K (T > 3K) for InAs (GaAs) 
QDs. These estimates for the relaxation rates due to 
one- and two-phonon processes are comparable to the 
ones recently measured in Refs. [1, [12], thus providing 
a reasonable explanation for these measurements. Note 
that, in contrast to the HH case, the relaxation time for 
electrons shows no deviation from the one-phonon time 
(or saturation) with decreasing i?-field 21 1 . 

Next, we provide explicit expressions of the relaxation 
rates for low and high temperature limits. The rates F^'^ 
can be written as 



Ft' it), 



(9) 



m=0 



where Si « 2'K{fi^eh\^(3'^ / K^ruhX^^Pccl), S2 « 
TTimtP^Et^/h'xy^cf), ds « TTimlP^Eyn^c^X^pl) 
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FIG. 1: The heavy-hole spin relaxation rate F for InAs QDs 
(GaAs QDs in the inset) as a function of magnetic field B for 
different temperatures T. The full curves represent the rate 
due to one- and two-phonon processes, i.e. F = X^f^i F'*' 
as defined in Eq. (|8} for difi'erent temperatures T, while the 
dotted lines present the one-phonon rate F'^'. 



= 5, r2 = 2, r3 = 0, and t — ksT/Eph with 
Eph = ^ci/\d. The functions F"''{t) depend on the 









F^''>{t) 


Fsit) 


t < 1 0.004 


0.015 


10»t" 10^*13 


5 • W'^t^'-^ 


10»ti5 


t > 1 0.08-!- 


0.03^ 


lO^t^ lO^t^ 


30*2 


0.3*2 


uJz 











TABLE I: The asymptotic values for -?;*'"■'(*). 

ratios t = ksT/Eph, d/X, and ci/ct. In Table |T]we list 
the asymptotic (scaling) expressions for Fj^"^\t) in low 
B-fields ujc.z ^ ^0 for low {t ^ 1) and high {t ^ 1) 
temperatures. We note that Fi'''\t) w Fi^\t) in both 
regimes, and F^^^'^'^^ = 0. 

Using Eq. ^ and Table [Jwe can write for the two- 
phonon rates, say, for InAs QDs 



F(2)=52 




0.5^) T<s:Eph 



+ 0.3^) T', T:>Eph 



f(3) = 5. 



109 ri5, T<^Eph 
0.3 T^, T:$>Eph. 



(10) 



From Eqs. (HO]) we find that for T < 2K and for B > 
0.5T the one-phonon processes dominate the relaxation 
rate F. On the other hand, for low B-fields (O.IT < 
B < IT) and finite temperatures (T > 2K) the two- 
phonon processes will give the main contribution to F, see 
Fig. [21 The main phonon processes could be identified 
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FIG. 2: The heavy-hole spin relaxation rate V in Eq. (|8]) 
for InAs QDs (GaAs QDs in the inset) as a function of tem- 
perature T for different B-field values. For finite B-field, F 
saturates at low temperatures due to one-phonon processes. 



experimentally by analyzing the temperature dependence 
of F, scaling as F ~ T for one-phonon processes and as 
F ^ fpj. two-phonon processes. Also, the saturation of 
F in vanishing B-field is a clear indication of two-phonon 
processes. Note that the strong enhancement of the two- 
phonon HH spin relaxation arises because (i) the rate is 
2nd order in SOI (whereas for electrons it is 4th order) 
and (ii) the effective mass for HHs is much larger than 
that for electrons. 

In order to compute F'^''^^ we took into account only 
the contribution from the deformation potential since 
this dominates the two-phonon relaxation for T/Eph > 
0.1 and iiJz,i^c ^ i^o- For the evaluation of F'^^ instead, 
we considered both the piezoelectric and deformation po- 
tential contributions, both of them being important for 
B and T considered here. Surprisingly, we found that 
the ZFR rate F^^) increases when decreasing the dot size 
as F*^'^) ~ A^^, while the other two rates decrease with 
decreasing the dot size as F(i) - and F^^) ^ A^. This 
behavior strongly differs from the electronic case where 
the ZFR mechanism is efficient for rather large dots [l^ . 

Interestingly, the present results do not change much 
if the B-field is tilted with respect to the QD plane. The 
g-factor for HHs is strongly anisotropic with <C gi_ so 
that one can neglect the in-plane Zeeman splitting. This 
implies performing the substitution uJc,z — > uJc,zcos9 in 
above results, with 9 being the angle between the B-field 
and the z-direction. This will lead to a reduction of the 
B-dependent rates (F^^'^))^ ^hile the ZFR (F^^)) being 
independent of B remains the same. 

In conclusion, we have shown that two-phonon pro- 
cesses give rise to a strong relaxation of the HH spin in a 
flat quantum dot. This time is predicted to be in the mil- 
lisecond range, comparable to the one measured in recent 
experiments on optical pumping of a HH spin in QDs Q . 



Though other sources of relaxation are not excluded, a 
careful scaling analysis of the measured relaxation time 
with the magnetic field and/or the temperature should al- 
low one to identify the two-phonon process as the leading 
relaxation mechanism for the heavy-hole spin localized in 
small QDs. 
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